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Abstract. The paper investigates cold molecules formation in the photoassociation 
of two cold atoms by a strong laser pulse applied at short interatomic distances, which 
lead to a molecular dynamics taking place in the light-induced (adiabatic) potentials. 
A two electronic states model in the cesium dimer is used to analyse the effects of this 
strong coupling regime and to show specific results: i) acceleration of the ground state 
population to the inner zone due to a non-impulsive regime of coupling at short and 
intermediate interatomic distances; ii) formation of cold molecules in strongly bound 
levels of the ground state, where the population at the end of the pulse is much bigger 
than the population photoassociated in bound levels of the excited state; iii) the final 
momentum distribution of the ground state wavepacket keeping the signatures of the 
maxima in the initial wavefunction continuum. It is shown that the topology of the 
light-induced potentials plays an important role in dynamics. 
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1. Introduction 

Photoassociation of cold atoms as a technique to form cold molecules has known 
new developments in the last years due to explorations using shaped laser pulses 
in order to control molecules formation: enhancement in cold molecules production 
and attainment of deeply bound vibrational states are both desired. The road from 
cold atoms photoassociation using continuous lasers [1] which keep the selectivity 
of transitions to the use of ultrashort pulses with broad bandwidth brought the 
challenges of a new physics. Prospective experiments using shaped femtoseconds 
laser pulses to photoassociate ultracold atoms [21 E] have shown the coherence of the 
process, but emphasized the difficulty to increase the number of created molecules. 
On the other hand, frequency-chirped light pulses in the nanosecond range were 
used to coherently control ultracold atomic Rb collisions showing the enhancement 
of short-range collisional flux [H Other experimental developments trying to 

establish coherent control techniques for cold molecules formation explored multiphoton 
photoassociative ionization in a Rb magneto-optical trap combining femtoseconds and 
continuous lasers [6] and optimal control of multiphoton ionization of Rb2 molecules 
using femtosecond laser pulses in a closed feedback loop [7] . Theoretical studies of pulsed 
photoassociation explored a variety of schemas to control cold molecular dynamics: 
with chirped pulses [8j El QUI E], pump-dump schemes to stabilize the cold molecules 
[T2| [T3l [HI [15], schemes using adiabatic passage [HI IT?] . 

The present paper prolonges previous works [HI [19] , the aim being to investigate 
theoretically the photoassociation of two cold atoms by strong laser pulses applied at 
small or intermediate interatomic distances. Such pulses have to be only "moderately 
strong" in order to avoid additional processes as ionization, and to act far from the 
atomic resonance (i.e. with a large red detuning, the colliding atoms being excited at 
small interatomic separations), to avoid the transfer of population to the continuum 
of the excited state. This regime of strong coupling and large detuning in cold atoms 
photoassociation can be used to address some specific interrogations: i) One interest 
is to explore if a strong pulse applied at small interatomic distances could be used to 
create strongly bound cold molecules, i.e. if such a pulse could accelerate efficiently the 
initial population, located at large interatomic distances in a cold collision, towards the 
inner region, ii) Secondly, a transition taking place at small or intermediate interatomic 
distances brings a regime necessarily "seeing" the nodal structure of the initial ground 
state continuum, and then we expect such traces in our results, iii) We are interested to 
explore the effects produced by a rather strong regime of coupling, generally not easy 
to be predicted. In [18] we have shown that for a strong coupling, characteristic times 
related to the adiabatic potentials become relevant in dynamics. The present work will 
explicitely explore the role played by the topology of the light-induced potentials in the 
dynamics and the photoassociation results. 

Our analysis is pursued on the example of the a 3 S+(6s,6s) — > l 9 (6s, 6^3/2) 
transition in CS2. The structure of the paper is the following: Section [2] describes the 



Cold molecules forr, 



iteraction 




R(a ) 



Figure 1. a 3 £+(6s, 6s) and l 9 (6s, 6^3/2) electronic potentials of CS2 (full lines), 
dressed with the photon of energy Tilul= Eq P3/2 — E 6s — HAl (HAl—I^O cm -1 ) 
and crossing in R c — 29.3 a , V c =Vi g (R c )=Vz(R c )= -143 cm' 1 . Dashed lines: the 
adiabatic potentials (V^ d , V^ d ) obtained from the diagonalization of the 2x2 potential 
matrix (fT3")l with the coupling Wl = 13.17 cm -1 . The energy origin is taken to be the 
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= of the l s (6s + 6^3/2) potential. 



time-dependent model and the time scales relevant in the photoassociation dynamics, 
as well as the initial wavefunction and the time evolution of the wavepackets during 
the pulse. In Section [3] we analyse the population transfer during the pulse and show 
the relevance of the light-induced potentials for the dynamics. In Section H] are shown 
and interpreted the results at the end of the pulse: formation of strongly bound cold 
molecules in ground and excited electronic states (Section 14.11) . acceleration of the 
ground state population to small interatomic distances (Section H~2l) . and the momentum 
structure of the final wavepacket in the ground state which reflects the maxima of the 
initial continuum wavefunction (Section 14.31) . An Appendix is connected to Section 14731 
Section [5] contains comments and conclusions. 

2. Simulation of the photoassociation dynamics during the pulse 

The photoassociation reaction is between two cold cesium atoms colliding in the 
ground state potential g = a 3 £+(6s,6s), at a temperature T = 0.11 mK, which are 
excited by a moderately strong laser pulse (with intensity I fa 43 MW/cm 2 ) to form a 
molecule in a superposition of vibrational levels {v e } of the excited electronic potential 
e = l 9 (6s, 6p 3 / 2 ). Only the s wave of the collision is considered. For a rotational 
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quantum number J = 0, the process can be schematized as: 

Cs(6s) + Cs{6s) + hu L -> Cs 2 (l g (6s + 6p 3/2 ); {v e }, J = 0) (1) 

The pulse is red-detuned with hAi = 140 cm -1 from the energy hu at = Eq Ps/2 — E 6s 
of the D2 atomic transition {hujL=huj a t- hA L ). The large detuning hA L determines a 
crossing of the field dressed diabatic potentials at the interatomic distance R c = 29.3 a , 
with V c =V lg (R c )=Vx(R c )=-U3 cm- 1 (figured]). 

The a 3 £+(6s, 6s) and l s (6s, 6p 3 / 2 ) electronic potentials used in the present 
calculation (figured!) are built from quantum chemistry [20] and asymptotic calculations 
[2~T| [22] and were described in a previous paper [18] . 




2.1. Time- dependent model 

The dynamics of the photoassociation process is simulated by solving numerically 
the time-dependent Schrodinger equation associated with the radial motion of the 
wavepackets ^i g (R, t) and *$?s(R, t) in the electronic channels l g and a 3 £+, coupled 
by an electric field with the amplitude £{t) = S f (t) cos u L t. The equation can be 
written as [T8] : 

\ = (2) 

*z(R,t) ) {) 

(T + V lg (R) W L f(t) \(y lg (R,t)\ 
\ W L f(t) f + Vx(R) ) \ *x(R,t) J 

Equation ([2D is obtained in the Born-Oppenheimer approximation for the diatomic 
molecule and using the rotating wave approximation with the frequency lvl/2tt. The 
potentials V\ {R) and V^(R) are the diabatic electronic potentials crossing in R c , 
represented in figure [TJ T is the kinetic energy operator and W^fit) the coupling 
between the two channels, with fit) the temporal envelope of the pulse. Wi = —^SqD^, 
where 8q = is the field amplitude (with I the laser intensity), e~2 the polarization, 
and the transition dipole moment between the ground and the excited molecular 
electronic states. We neglect the R-dependence of the transition dipole moment, using 
the asymptotic value D e % deduced from standard long-range calculations for a linear 
polarisation vector el. This approximation remains good for the present calculation, 
as for distances around the crossing R c = 29.3 a , the dipole moment is closed to its 
asymptotic value, D e g \[R c ) rs 0.9 D e g \ , and decreases slowly for smaller distances [23J. 
For a pulse intensity I ~ 43 MW/cm 2 , the coupling becomes Wl=13.17 cm -1 , inducing 
a significant avoided crossing, as it can be seen in figure [Tj, where the light-induced or 
adiabatic potentials V^ d , V£ d are represented with dashed lines. 

The numerical calculations were made for a rectangular pulse with a duration of ~ 
300 ps, whose envelope f(t) is represented in figure El 

The Schrodinger equation ([2]) is solved by propagating in time an initial 

wavefunction j 5* , ) on a spatial grid with the length L R = 760 a . The time 
\ *s(.R,0) ) 
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Figure 2. Temporal envelope f(t) of the photoassociating pulse. 

propagation uses the Chebychev expansion of the evolution operator [2H EH] and the 
Mapped Sine Grid (MSG) method [TUJ [27] to represent the radial dependence of the 
wavepackets. 

The results extracted from the dynamics are: 
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• the evolution of the wavepackets during the pulse, for the two channels g = a 3 S 
e = l g , in the position representation, "$?-£,ig(R, t), and momentum representation, 
defined by the Fourier transforms ^^^ g (p,t); 

• the evolution of the population in each electronic state during the pulse. At a given 
instant t, the population in one of the electronic states g, e is calculated on the 
spatial grid extending from R m i n to L R , as: 

p g , e (t)= [ Lr \^ g AR',t)\ 2 dR' (3) 

The spatial grid is chosen such as at every instant t the total population is 
normalized at 1 on the grid (-Ps(i) + -Pi (t) = 1)- At t=0 the population is entirely 
in the ground state (-Ps(O) = 1). 

2.2. Time scales related to the laser coupling and vibrational movement 

The time scales relevant for the dynamics are related to the laser coupling and to the 
vibrational movements in the electronic potentials. The spontaneous emission from the 
excited state is neglected, as the time evolution of hundreds picoseconds studied here is 
short compared with the spontaneous emission time of about 30 ns. 

We begin by defining the characteristic times connected to the laser coupling. 

A local time- dependent Rabi period can be associated with the laser coupling Wi,f{t) 
between the two electronic states [9]: 

TjudR, t) = hn (4) 

'(W L /(t)) 2 + A2(J2) 



where 2 A (R) = \Vi g (R)—Vs(R)\ is the local detuning (see figured]). Such a characteristic 
time is relevant if the impulsive approximation remains valid on the whole duration of 
the pulse, i.e. if the relative motion of the two nuclei can be considered as frozen during 
the pulse duration. 
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For the rectangular pulse studied here the coupling remains constant in the time 
interval (100 ps, 300 ps), so we can refer in the analysis at a local Rabi period associated 
to the constant coupling Wl- 

Trou(R) = 1 ^ (5) 



This local Rabi period has its maximum at the potentials crossing (T^ abi (R c )=l .27 ps), 
diminishing with the increasing of the local detuning (for example T^ abi (89 ao)=0.24 
ps). 

One can also associate a Rabi period with the beating induced by the coupling 
Wl between two specific vibrational states, one belonging to the excited electronic 
state, and the other to the ground state. Indeed, in equation ([2j) the wavepackets 
can be developed as superpositions of vibrational wavefunctions {\Xv e ,v g {R) >} with 
eigenenergies E VejV , corresponding to each electronic Hamiltonian H eg = T + V e g 
(H e ,g|xu e ,i! g >= E VetVg \xv e ,v g >)• 

*i g (R,t) = E c ^)exp(-^ e t)x, e 0R), (6) 

MR,t) = Y, c v a ( t ) e M-TE Vg t) Xvg {R), (7) 

Vg n 

Supposing only two vibrational states, c Ve (t) exp(—jrE Ve t)x Ve (.R) and c v (t) exp(—±E v t) 
Xv g (R), respectively, with \c Ve {t)\ 2 + \c Vg (t)\ 2 = 1, (^(O)] 2 = 1, (^(O)] 2 = 0, which are 
coupled by Wl in equation (T5]), one obtains for the oscillating population in the excited 
state: 

= ttt^ d sm {Q VetVg t), (8) 

{ill iv e ,Vg) 

^ Ve , Vg = ^\W L < X v e \Xv g > I 2 + i(E Ve - E Vg )/2Y. (9) 

The corresponding Rabi period will depend on the overlap < X v e \ X v g > of the vibrational 
functions: 

T L = (10) 

Ve,Vg Q V iW ; 

iL V e ,Vg 

The time scale associated with the vibrational motion of a vibrational level v with 
binding energy E v in an electronic potential is: 

Tf= , 27Tk , . (11) 

To estimate the coupling influence on the dynamics, the time scales (JSJ) and (flOl) 
associated with the laser coupling shall be compared with the vibrational periods (TTTT) 
of levels in the ground and excited potentials. In the case treated here, the Rabi periods 
associated with the laser coupling are of the order of picosecond, being much smaller 
than the characteristic vibrational periods implied in the problem (tens or hundred ps), 
which indicates of strong coupling [IB] . 
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Figure 3. Initial state of the photoassociation process of energy Eq: stationary 
continuum state in the a 3 H^ potential, calculated and normalized to 1 in a box of 
length Lr w 760 ao, corresponding to the temperature E o /kB=0.ll mK. (a) Position 
representation \^^(R, 0)|. (b) and (c) Momentum representation \^Y,{k,0)\ shown in 
the domains of "large momenta" (\k\ < 20 a.u.) and "small momenta" (|fc| < 0.2 a.u.), 
respectively. 

2.3. Initial state: spatial and momentum representations 

The representation of the initial state in a wavepacket treatment of the cold atoms 
photoassociation is discussed in [10]. For a low temperature collision (T =0.11 mK) and 
an excitation process taking place at short distances (R c = 29.3 ao), the initial state of 
the photoassociation process has to be represented using stationary collision states. We 
shall simulate the photoassociation dynamics using as initial state a continuum state 
of the ground electronic channel, having the energy E = k^T . The results obtained 
can be used to estimate a photoassociation probability for an ensemble of cold atoms in 
thermal equilibrium at temperature T (see Section [4]). 

In our numerical method, the initial continuum state is calculated as one of 
the eigenstates of the ground electronic state Hamiltonian, through the Sine Grid 
Representation [27], in a box of length Lr. The method introduces a discretization 
of the continuum which supply continuum states having a node at the boundary of 
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Figure 4. Time evolution of the wavepackets in the position representation: \^>s(R, t)\ 
(left column) and \ 1 ^i g (R, t)\ (right column). The vertical dashed line in the left 
column marks the part of the a 3 S+ wavepacket which is strongly accelerated inside 
the potential during the pulse. 



the box (as the sine basis functions used in representation). Then it is possible to 
have a continuum delocalized state as initial state in the wavepackets simulation of 
the photoassociation dynamics. The MSG method allows the use of large spatial grids 
on which the wavepackets dynamics in the range of distances relevant for the problem 
can be followed for long propagation times without being influenced by the external 
boundary of the box [TO] . 

Here the initial state ^^(R, 0) (see the (a) panel in figure [3]) is chosen to be a 
continuum wavefunction of the a 3 £+(6s,6s) potential, of energy Eq = 7.6908 x cm" 1 
corresponding to a temperature T = E /k B =0.11 mK. This wavefunction is calculated 
through the MSG method [271 HO] in the box of length L R w 760 a , with a node at the 
end of the grid, and normalized to 1 on the grid. To obtain the normalization per unit 
energy the populations have to be multiplied by the density of states in the box at Eq, 
[{dE / dn^EoY 1 [10J. The energy resolution for neighbouring eigenstates in the box at 
E is 5E\ E=Eo = {dE/dn)\ Eo = 0.632 x 10" 4 cm" 1 , corresponding to [5E\ E=Eo ]/k B =0.09 
mK. 

Figures [3]^b) and (c) show the amplitude of the initial wavefunction in the 
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Figure 5. Time evolution of the wavepackets in the momentum representation 
*s(fc,t)| (left column) and \^i g (k,t)\ (right column). 



momentum representation, |\l/s(fc,0)|. The wavefunction amplitude is mainly localized 
in the domain of "small momenta", \k\ < 0.06 a.u. (figure E^c)), which corresponds to 
the large distances domain (R > 100 ao) in the position representation. A picture from 
the domain of "large momenta" (\k\ < 20 a.u.), where the wavefunction amplitude is 
much smaller, is displayed in figure [3]^b). This domain of momenta corresponds to the 
domain of "small distances" in figure [3](a). We are interested to observe the changes 
appearing in both domains during the time evolution. The negative k values mean 
momenta oriented to the inner wall of an electronic potential, and the positive k values 
momenta oriented to large distances. 



2.4- Wavepackets evolution during the pulse: spatial and momentum representations 

The evolution of the wavepackets during the pulse is illustrated in figures H] and El 
showing the a 3 S+ and l g wavepackets in the spatial and momentum representations, 
respectively. The dynamics will be analysed in order to understand the vibrational 
movements inside each channel and the exchange of population between the electronic 
channels (figure [6]). 

In the l g potential, the excited wavepacket extends on the whole spatial grid, 
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showing that bound and continuum levels are populated during the pulse. The 
vibrational movement of the population occupying bound states with outer turning 
points in the crossing region (R c = 29.3 do) can be well observed in the right columns 
of figures H] and [51 The large amplitude of the wavepacket in the zone of big momenta 
(t=170, 250, 295 ps in the right column in figure HI) is equivalent with a strong presence 
of population at R < R c at the same moments. After the pulse, only bound levels remain 
populated in l g (figure [8] and section HJ). 

In the a 3 S^ potential, the wavepacket moves progressively to the inner region 
(left column of figure HJ). The vertical line in the same figure marks the separation 
in two spatial domains, which are differently affected by the pulse: at small and 
intermediate distances R < 70 a , the wavepacket is accelerated and deformed, but 
at large interatomic separations R > 80 ao the action of the pulse can be considered as 
impulsive. 

The wavepackets dynamics in the momentum representation (figure [5]) makes visible 
an unexpected feature which appears from t=150 ps in the time evolution of the ground 
state wavepacket ^^(k,t) at k « 12.2 a.u., and whose intensity increases until the end 
of the pulse. This result will be analysed in section 14.31 

3. Analysis of the population transfer during the pulse 

We shall analyse the time evolution of the population transferred by the pulse in 
the excited state, Pi g (t), calculated with formula ([3]) and displayed in figure M (the 
population in the ground state has a complementary evolution, as -Ps(t) + Pi g (t) = 1, 
with -Ps(O) = 1). The figure shows that from the large amount of population Pi g {t) 
transferred during the pulse, which is of the order of 8.7 x 10~ 3 , only 7.77 x 10~ 5 
remains at the end of the pulse. As it was discussed in the previous section, the large 
population transfer in the excited state during the pulse is related to the occupation 
of continuum states at large distances, in levels which do not rest populated after the 
pulse [9j [ID] • On the contrary, what it is interesting for us is the population in bound 
vibrational levels, counting as cold molecules formation. 

The vibrational levels of l g and a 3 £„ predominantly populated by the pulse were 
identified by calculating the probabilities P Ve ,v g (t) = \ < ^ e , g (R,t)\xv e ,v g (R) > | 2 that 
a certain vibrational level of the excited or ground electronic state (v e or v g ) to be 
populated at an instant t. Two kinds of bound vibrational levels are mainly populated 
during the pulse in the l g state: 

• levels around the crossing of the diabatic potentials: t> e =141-144, with vibrational 
energies E Ve =-144.5, -141.39, -138.29, -135.24 cm -1 , lying in a domain of about 
10 cm -1 containing the crossing energy V c = —143 cm -1 , and with outer turning 
points of the vibrational wavefunctions lying around R c = 29.3 ao- These levels 
have vibrational periods of about 11 ps (see tabled]). 

• levels with v e > 241 (E Ve=2 4i = —8.818 cm -1 ), having vibrational wavefunctions 
lying at much larger distances (R > 60 a ). From these levels populated at large 
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Figure 6. Time evolution of the population in the excited state, P\ g {t), in the 
photoassociation with the pulse Wi,f(t) (whose envelope f(t) is represented in figure^]). 
The inset shows the details of the time evolution during the constant coupling Wl 
(between 100 ps and 295 ps). 



distances, only the levels v e =244 up to 248 remain notably populated after the 
pulse. Their vibrational periods are in the range 116 up to 130 ps. 

The inset of figure [6] shows the evolution of P\ (t) in the time interval of constant 
coupling Wl (between 100 ps and 295 ps). The oscillating features of Pi g {t) during 
this time interval come from the exchange of population between vibrational levels of 
l g and a 3 E+ located around the crossing (t> e =141 up to 144 in the excited state, and 
f 5 =44 up to 49 in the ground state, see tabled]), without major contribution from the 
population transferred at large distances. Indeed, as it is marked on the inset of figure [6] 
the energy domain covered by these oscillations is about 15 x 10 -5 , comparable with 
the final population Pi (t = 395ps) = 7.77 x 10~ 5 in bound levels. Also, the variations 
of the total population Pi g (t) between instants as t=130, 150, 170 ps, etc. are very 
closed to the variations of the population in levels around the crossing, given by the 
sum P lg ~ (t)=X^ e =i4i Pv e (t)- Comparatively, the probability for the population of 
the levels v e =244-248, P? g 4 ~ M8 (t)=J2 2 v % 44 P Ve (t), shows much smaller variations. 

The quantum beats appearing in figure [6] are related to the characteristic times 
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Table 1. Vibrational levels v e of the excited state (with energies E Ve and vibrational 
periods T" b ) and v g of the ground state (with energies E v and vibrational periods 
Ty lb ), populated during the pulse and giving the quantum beats shown in the inset 
of figure [6] {v e } are l g vibrational levels located around the crossing of the diabatic 
potentials and {v g } are vibrational levels in the o 3 S+ potential giving the biggest 
overlaps < Xv g \Xv a > with every v e level. The last two columns show the energy 
differences \E Vg —E Ve | and the characteristic times Ue , calculated with formula (I10| . 
The energy origin is the dissociation limit Eq s+ q P3/2 of the l ff (6s, 6^3/2) potential. 
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of the dynamics: the vibrational movement in each potential well and the beating 
between the two coupled wavepackets, which are superpositions of vibrational functions 
corresponding to each electronic potential, as in (jBJ) and (JTj). Table [H gives a list of 
energies E v , vibrational periods T mb , and characteristic beating times Ve between 
levels in the ground and excited electronic states which contribute significantly in the 
exchange of populations between the two coupled channels; it also contains the overlaps 
< Xv a \Xve > with a 3 E+ levels. Only the v g levels having the biggest overlaps with a given 
v e level are shown: these are the levels t> 9 =44 up to 49, whose vibrational wavefunctions 
have the outer turning points at distances R < 35 ao in the a 3 S+ potential, and which 
are strongly populated by the pulse. 

The population transfer between the two electronic channels is regulated by two 
time scales: a longer one, related to the vibrational movements inside the potential wells, 
and a shorter one, related to the Rabi coupling. The longer scale reflects the influence 
of the vibrational movements on the exchange of population between the two channels: 
the exchange is maximal when the amplitudes of the two wavepackets have a significant 
overlap, which, in the case of wavepackets vibrating in two different potential wells, 
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arrives when both wavepackets have important localization probabilities in the crossing 
region [18J . When one of the wavepackets vibrates inside its potential, the transfer is 
generally diminished. The levels in the ground state have vibrational periods (> 40 ps) 
much longer than the levels in the excited state (about 11 ps). This explains the period 
of 40 ps appearing in the P\ (t) oscillations, which coincides with the vibrational period 
of the level t> 9 =44 in the a 3 S+ ground state, whose energy -143.27 cm -1 is very close to 
the crossing energy V c . On a much shorter scale, the transfer is guided by the strong 
laser coupling, which couples differently the implied levels. Table [1] shows characteristic 
times Ve of beating (calculated with formula [TOl) varying from 2.5 to 12.7 ps for 
coupled pairs of vibrational levels in the ground and excited states. Comparing these 
times with the vibrational periods of the concerned levels, the strength of the coupling 
appears as varying very much among pairs of levels. The times scales given by Ve 
are indicative for the short Rabi times appearing during the dynamics, for example the 
periods of 3 up to 4 ps of the fast oscillations in Pi (t). 

During the pulse the population accumulates in the l g state, such as in figure [6] 
appear not only the population beatings between the two channels, but also traces of 
the vibrational dynamics inside the excited potential: around t=250 ps, the two peaks 
of Pi g {t) are separated by a time interval of ~ 11 ps, which is the vibrational period of 
the v e = 142, 143 levels located around the crossing: T™^ 142 143 11 ps. 

3.1. The light-induced (adiabatic) potentials 

The mechanism of the population transfer during the pulse is enlightened if one 
analyses the light-induced (adiabatic) potentials. In figure [TJ are represented both the 
diabatic potentials Vi and the adiabatic ones, V^ d and V^ d , obtained from the 

diagonalization of the 2x2 potential matrix V el with constant coupling Wl as non- 
diagonal term: 



(Ve,+w l )u= y > ", m (12) 




(V. 1+ W t )U S =(^ (fl) V £ (R) ) (13) 

Expressions (1121) and ( {TBI illustrate the diabatic and adiabatic representations, 
respectively, of the V e i + W L operator. The coupling produces a strong deformation of 
the diabatic potentials Vi g , Vy, around the crossing (figured]), but its influence goes far 
beyond the crossing region, and this can be clearly seen if one compares the vibrational 
levels of the diabatic potentials V\ , Ve with those of the adiabatic ones V^ d , V^ d . Such a 
comparison can be made using the series of rotational constants B v = (xv\h 2 / (2[iR 2 )\xv) 
characterizing every potential. The energies E v and rotational constants B v were 
computed by solving numerically the corresponding stationnary Schrodinger equation 
(T + V)|x„ >= E v \xv > using the Mapped Fourier Grid Hamiltonian (MFGH) method 
[26J. In figure [7] are represented the rotational constants B v as functions of the 
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Figure 7. Rotational constants B v = (x v \h /(2fj,R 2 )\x v ) a s functions of the 
vibrational energies E v , for vibrational levels of the l g and a 3 S^ diabatic potentials 
(full and empty circles, respectively), and of the adiabatic potentials V£ d , V^ d (stars). 
The vertical line at V c = —143 cm" 1 indicates the crossing region. 



vibrational energies E v for a 3 £+ and l g diabatic potentials, as well as for the adiabatic 
potentials V^ d , V^ d . We use the B V (E V ) functions to observe how the energies of the 
vibrational levels in these potentials are distributed in a domain lying between -170 
and -128 cm -1 around the crossing energy V c =-143 cm -1 . These results show that the 
influence of the coupling is strongly felt in a large energy domain of several tens of cm -1 
around the crossing. 

Also marked in figure [7] are the vibrational levels populated during the pulse and 
relevant to the dynamics, together with their vibrational periods. We shall focus on the 
levels located around the crossing region. The levels v g = 44, 45 in the a 3 £„ potential 
have vibrational periods T™ b = 40, 49 ps. In the same energy region, the levels of the 
Vad potential have bigger vibrational periods T™ a b d = 43, 58 ps. Also, if the vibrational 
period of the l g levels located around the crossing, v e = 141 up to 144, is T™ b ~ 11 
ps, in the adiabatic potential V^ d the levels belonging to the same energy domain have 
a vibrational period twice bigger: T^ d w 23 ps. Then, as it could be expected from 
the shape of the V^ d and V^ d (figure [TJ, the vibration in the adiabatic potentials is 
slowed down in the crossing region, influencing the population transfer between the two 
channels. The slowing down of the vibrational movement in the excited state lead to 
longer periods for the population exchange between the two channels and constitutes a 
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Figure 8. l g wavcfunction at the end of the pulse. 



mechanism for increasing progressively the population in the excited state. 

The analysis of the adiabatic potentials can be used to extract characteristic times 
for the population exchange between the two channels. In an analogy with a two 
coupled states system, the Rabi beatings in the exchange of population is related to 
the Bohr frequency of the coupled system, which here can be found from the typical 
frequencies \E e ad — E® d \/2h, connected with the new energies E^ d , E 9 ad of the levels in 
the coupled system. In figure [7] we show that for levels located in the crossing region 
such a characteristic time T osc = 2irh/\E a \ d — E 9 ad \ is about 45 ps, close to the period of 
~ 40 ps of the large oscillations in figure [6j 

4. Results at the end of the pulse 

Figures [S] and M show the wavepackets l g (in position representation) and a 3 S+ (in 
position and momentum representations) at the end of the pulse (t=395 ps). The main 
characteristics of the results will be discussed in the following. 

4-1. Formation of strongly bound cold molecules in a 3 S+ and l g electronic states 

The strong coupling between the ground and excited states creates an interesting result 
at the end of the pulse: the population in bound levels of the ground state (mainly 
six levels, v g =47 up to 52), P^ 7 ~ 52 (t=395 ps)=2.83 x 1(T 4 , is much bigger than the 
population photoassociated in the excited state, Pi ff (t=395 ps)=0.78 x 10~ 4 . Some of 
the vibrational levels populated in a 3 S+, for example v g =47,48,49, have wavefunctions 
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Figure 9. a 3 S+ wavepacket at the end of the pulse, t=395 ps. (a) Position 
representation \^>^(R,t = 395ps)|. (b) and (c) Momentum representation \^j^(k,t = 
395ps) \ in the domain of "large momenta" (|fc| < 20 a.u.) and "small momenta" (|fc| < 
0.2 a.u.), respectively. 



localized at distances R < 35 a , being then strongly bounded. 

The final population in the 1 5 excited state (figure E]) is entirely in bound vibrational 
states. A superposition of two kinds of vibrational levels is created, showing two 
mechanisms in the population transfer: at resonance, where mainly one level rests 
populated, v e =143 (with outer turning point at ~ 30 ao), whose population represents 
82% P lg (t=395 ps), and off-resonance where several vibrational levels with outer 
turning points around R ~ 70 a (t> e =244 up to 248, representing 15% Pi 9 (t=395 
ps)) are populated due to the strong coupling catching the large amplitude of the initial 
wavefunction between 45 and 65 ao (figure [T0l) . This last kind of transfer creates a 
hole in the ground state in this domain of interatomic distances, as it is indicated in 
figures My&) and [TQ1 Such a result reflects the specificity of the present photoassociation 
conditions of strong field and large detuning. 

The results obtained for a total population normalized at 1 on the grid allow an 
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estimation of the averaged probability corresponding to a thermal distribution [TU] at 
the temperature T = E Q /kB=0.11 mK, as: 

P^m^F^iEo^E,)- 1 ^- (14) 

where Pi 9 (£ , )=Pi fl (t=395 ps)=0.78 x 10~ 4 is the l g probability obtained in the 
present calculation with an initial state of energy E , [(d£/dn)| So ] _1 is the density 
of states in the box of length Lr at Eq, ks is the Boltzmann constant, and Z is the 
partition function for a gas composed of non-interacting pairs of atoms in a volume V 
(with /j, the reduced mass of the diatom): Z = (2nfik B T) 3 / 2 V/h 3 . We then obtain 
ZP lg (0.11 mK) = 0.95 x 10~ 4 and ZP^ 3=47 ~ 52 (0.11 mK) = 3.44 x 10" 4 . For a 
number of N atoms in a volume V, and taking into account the spin degeneracy of the 
Cs(6 2 S) atomic state and of the initial electronic state, the total number of molecules 
photoassociated in l g per pump pulse is [TU]: N mol) i g = ^y-P lff (T)|. For a volume 
\^=10 -3 cm 3 and a density of atoms N/V=10 u atoms/cm 3 , the number of molecules 
obtained at the end of the pulse are: N mo i t i g = 0.3 xlO -2 and N mo i^,v „=47-52 = 1-lx 10~ 2 . 

4-2. Acceleration of the a 3 S+ wavepacket to small interatomic distances 

During the pulse, the a 3 E+ slow packet is accelerated towards small interatomic 
distances R. The left column of figureHJshows the time evolution of the ct 3 £+ wavepacket, 
which advances progressively to the inner zone. Especially the part of the wavepacket 
occupying distances R < 70 a feels the acceleration to the crossing zone, where the 
diabatic potential a 3 S+ in R~ 6 becomes the adiabatic one V^ d decreasing in R~ 3 (see 
figure [lj. On the contrary, the maximum of the wavefunction located at R m 89 ao 
does not move during the pulse, but begin to be accelerated after the pulse. Similar 
observations can be made on figure [9] (a), showing the a 3 S+ R- wavepacket at the end 
of the pulse (t=395 ps): the changes in the wavepacket amplitude are at distances R 
< 89 ao- The momentum representation of the a 3 S^ wavepacket at the same instant 
t=395 ps, in figure [9] (c), shows that, compared with the initial symmetric distribution 
(figure [3] c) , the part corresponding to k<0 has moved to bigger \k\ values, which also 
indicates the gain of kinetic energy in the electronic potential in the movement to small 
distances. 

The creation of a "hole" in the ground state wavepacket at the end of the pulse 
(discussed previously and shown in figure [9^)) is also a factor leading to a compression 
of population at short range, which after the pulse acts to increase the density of atom 
pairs at short distances [TT] . 

4-3. Kinetic energy "gains" in the final a 3 S+ wavepacket as signatures of the maxima 
in the initial wavefunction continuum 

In this section we shall analyse some peculiar features appearing in the ground state 
wavepacket \^?z(k, t) \ during the time evolution and at the end of the pulse (figure [9b)). 
The time evolution in the momentum representation (left column of figure [5]), shows 
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Figure 10. Details of the initial wavefunction on the ground state, \^>^(R, t — 0)|. The 
arrows show local maxima of the wavefunction. The encircled ones give contributions to 
the off-resonance cycling of population between the two channels, bringing momentum 
to the final ground state wavepacket. On the contrary, ground state population from 
the local maxima at R=47 ao and 59 ao is tranferred to the excited state l g (figure [8]), 
creating a hole in the final a 3 E+ wavepacket (figure [9^i) . 

that, from t=150 ps, a line strikingly appears in the wavepacket amplitude \^Y,(k,t)\, 
at k ~ 12.2 a.u. The kinetic energy associated with this momentum is h 2 k 2 /2/i= 134.8 
cm _1 = \Vi g (R) — VE(-R)|_R=89a 5 corresponding to the local difference in energy between 
Vi (R) and Vs(R) potentials at R = 89 a . The amplitude of the initial wavefunction 
|\I/e(-R, 0)| has a maximum at this distance (figure [3] a), which does not move during 
the time evolution, but on which small oscillations begin to be superposed (see the 
left column of figure H]). The period of these oscillations is T R = 0.5 a = 2n/k, 
corresponding to a plane wave exp(ifci?) with k = 12.2 a.u. The momentum width of 
this fc-feature in the i)| wavepacket is Ak « 0.15 a.u., which is consistent with 

the spread in distance AR > 3.3 a on \^/^{R, t)\ around R ~ 89 ao- We interpret 
this feature as corresponding to the kinetic energy of the population cycled back in the 
ground state from the excited state in the fast exchange of populations taking place 
around R ~ 89 ao due to the strong coupling between the two electronic channels, 
with Tft abi (89 a )=0.24 ps. At the end of the pulse, the line with k w 12.2 a.u. is 
accompanied in the \^s(k,t = 395 ps)\ wavepacket by its negative value k = —12.2 a.u. 
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(see figure [9b), which has to be a technical artefact coming from the propagation of 
this large positive momentum component to large interatomic R distances, followed by 
reflection at the end of the grid. 

In fact, as it is shown in the Appendix, it seems that the momentum k first appears 
associated with an ingoing plane wave travelling to small distances (k < 0), but it is fast 
reflected by the inner wall of the potential and transformed in a positive k travelling 
to large distances, which is easy to be observed in the wavepackets evolution. We have 
simulated the dynamics in the same conditions, but taking as initial wavepacket in the 
ground state a gaussian centered at R = 89 a , and indeed we observed a peak at 
k = —12.2 a.u. appearing early in the time evolution. 

k = 12.2 a.u. is not the only value of momentum for which a feature appears in 
the i)| wavepacket. At various instants of the time evolution, other lines can 

be observed, at smaller values of k. As these lines are embedded in the wavepacket, 
it becomes easier to distinguish them for larger momenta where the wavefunction 
amplitude is smaller. In figure [9b, showing the a 3 £„ wavepacket amplitude in the 
domain of large momenta, at the end of the pulse, other lines can be observed at the 
values: k=12.2, 10.4, 8.6 and 6 a.u. They correspond to kinetic energies « 2 fc 2 /2/i=134.8, 
98, 68.6 and 32.6 cm -1 . If one associates these kinetic energies to local differences 
AV(R) between the electronic potentials, then, according to the reasoning just exposed 
they correspond to local transitions from the excited to the ground state taking place 
at distances R ~ 89, 43.9, 36.8 and 32 ao, where maxima of the initial wavefunction 
are located, as it appears in figure [101 Other local maxima of the initial wavefunction, 
at R=47 ao and 59 ao, do not have a correspondent k- value in the final momentum 
distribution of the ground state wavepacket, but in this domain of distances the 
population is not cycled back to the ground state, remaining transferred to the levels 
f e =244 up to 248 of the excited state (see figured]). As a consequence, as we showed 
already, a hole is created in the ground state in this spatial domain (figure [9k) • 

We also have to mention that for longer propagation times one can see appearing 
in the \^z(k,t)\ wavepacket peaks corresponding to the negative k values of these other 
smaller momenta: -10.4, -8.6 and -6 a.u. 

In the Appendix we use the impulsive approximation in the limit A(i?) ^> Wl 
to show the emergence of such momenta during the time evolution. Indeed, as we 
emphasized, for the maximum of the ground wavefunction located at R = 89 ao, the 
impulsive approximation rests valid during the whole pulse (see figure Eh), and for other 
maxima the impulsive approximation could be applied on smaller durations of the pulse. 

The significant fact is that these "large momenta" appearing at the end of the pulse 
in the ground state wavepacket are signatures of the maxima in the initial wavefunction 
continuum. This is a specific effect of the strong photoassociation pulse applied at small 
distances which reveals the structure of the initial ground state. 
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5. Comments and Conclusions 

We have analysed the dynamics in the photoassociation of a pair of cold atoms 
by a strong laser pulse (I ~ 43 MW/cm 2 ) applied at short interatomic distances 
(R c = 29.3 a ) for a cold collision. The numerical calculations were made for the 
a 3 E+(6s + 6s) — > l«y(6s + 6p 3 / 2 ) transition in Cs 2 , at a temperature T = 0.11 mK. 
The large detuning (hA L = 140 cm -1 ) and the intensity of the pulse were chosen to 
correspond to a specific regime imposed by the limit of an asymptotic detuning much 
bigger than the maximum of the coupling, hA^ = 140 cm -1 ^> Wl = 13.17 cm -1 , in 
order to avoid the transfer of population to the continuum of the excited state at the 
end of the pulse. 

The specificity of this regime comes from two sides: a) the large detuning which 
locates the resonance condition at small or intermediate interatomic distances, making 
"visible" the nodal structure of the initial continuum wavefunction; b) the strong 
coupling between the two electronic channels, which acts not only on vibrational levels 
around the crossing of the dressed electronic potentials, but also induces off-resonance 
cyclings of population between the coupled electronic states. 

We have chosen not only a quite strong pulse, but also a quite long one (a 
rectangular pulse of about 300 ps), in order to obtain a better understanding of the 
time evolution in the presence of the pulse and to see how efficient for the cold molecule 
formation is the acceleration of the population from large to small interatomic distances. 

In this strong regime of coupling, the photoassociation dynamics during the pulse 
takes place in the light-induced (adiabatic) potentials, whose topology influence the 
exchange of population between the coupled channels. In our example, the shapes of 
the adiabatic potentials lead to acceleration of the ground state population to the inner 
region (at short distances the diabatic potential a 3 £+ in R~ 6 becomes the adiabatic one 
V^ d decreasing in R~ 3 ), and also to the slowing down of the vibrational movement in 
the crossing region, which is a mechanism for increasing progressively the population in 
the excited state. 

The main characteristics of the results at the end of the pulse are the following: 

(i) It appears that such a pump scheme allows for the production of ground state 
molecules through a single laser pulse. Indeed, strongly bound cold molecules are formed 
in a 3 S+ and l g electronic states, the population transferred in bound levels of the ground 
state being much bigger than the population photoassociated in bound levels of the 
excited state. Some of the vibrational levels populated in a 3 S+ have wavefunctions 
localized at distances R < 35 ao. The final population in the l g excited state is entirely 
in bound vibrational levels, populated both at resonance (v e = 143, with outer turning 
point at pa 30 ao) and off-resonance, where several vibrational levels with outer turning 
points around R pa 70 ao are populated due to the strong coupling catching the large 
amplitude of the initial wavefunction between 45 and 65 a . This last kind of transfer 
creates a hole in the ground state in this domain of interatomic distances. 

(ii) During the pulse the population in the ground state is globally accelerated 
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towards small interatomic distances. The creation of a "hole" in the ground state 
wavepacket at the end of the pulse is also a factor leading to a compression of population 
at short distances [TT] . 

(iii) At the end of the pulse, the momentum distribution of the ground state 
wavepacket keeps the traces of the initial continuum maxima. This is an effect of 
the strong coupling leading to off-resonance cycling of population between the two 
channels and bringing kynetic energy in the ground state. The cycling of population is 
particularly important at those interatomic distances where the maxima of the initial 
continuum are located. 

An important question is if the regime explored here can be mantained if, for 
example, one increases the detuning HAl, in order to form cold molecules in lower 
vibrational levels of the ground state. An insight about the results which can appear 
by increasing the detuning can be extracted in the impulsive limit, which normally rests 
valid at large interatomic separations. In the impulsive limit, the population in the 
excited state at the end of the pulse can be approximated as [21] : 

W 2 

|Vl/e(jM|2 ^ Wi + A\R) S ™ 2 (™)\%(R,0)\ 2 , (15) 

where hQ(R) = ^/Wl + A 2 (R) (for R^oo, 2A(R) -> HA L ). Equation (HHJ) shows that 
the transfer from the ground to the excited state is favoured: a) at large interatomic 
distances R, because \^ g (R, 0)| 2 is the amplitude of an initial continuum wavefunction 
of low energy in the ground state; b) at resonance (A(R) = 0); c) for a strong coupling 
Wl, when, depending on the ratio Wl/JiAl and the pulse duration, various ranges of 
distances can be populated. It has to be observed that, increasing the detuning A(R), 
one increases fl, so the oscillations in sin 2 (fit) with the period tt/Q will be very fast 
and produce an unstable regime of transfer, in which the result at the end of the pulse 
cannot be predicted. In this regime, the large distances and the continuum of the excited 
state can easily rest populated after the end of the pulse. The only manner to overcome 
these difficulties is to increase the coupling, which can broke the impulsive limit and 
bring population to short distances. Then, the ratio Wl/TlAl is indeed a significant 
parameter for the regime explored here, which can be maintained for a larger detuning 
only if the pulse intensity is also increased. 

The acceleration of the population to the inner region and the efficiency in forming 
strongly bound cold molecules is related here to a non-impulsive regime of coupling 
brought at R < 70 ao by the strength and the time duration of the pulse. 

These results have to be completed with a further investigation of the time evolution 
subsequent to the pulse. Also the model can be enriched by considering couplings with 
other electronic surfaces. Our work is a tentative to show that the photoassociation 
of cold atoms at small/intermediate distances with an intense laser pulse offers the 
possibility to explore cold molecules formation using the specific topologies of the light- 
induced potentials. Favorable conditions can be found taking into account the variety 
of electronic transitions in alkali dimers which can be controlled by the parameters of 
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Figure Al. The momentum representation \^g A (k,t — 300 ps)\ of the wavefunction 
in the ground state ^g A (R,t), calculated in the impulsive approximation with (|A.1|) . 
for t=300 ps. The initial wavefunction <3> ff (i?,0) has maxima at R = 33, 36.8, 41, 
47, 59, 89 ao (see figure [TU]), which, according to our reasoning, are supposed to give 
peaks in the momentum distribution ^/g A (k,t) of the ground state, at k ss 6, 8.6, 
10.4, 10.8, 11.6, 12.2 a.u. The figure shows that indeed the momentum distribution 
calculated in the impulsive approximation displays such peaks, at a time t close to 
the pulse duration. It can be noted that the amplitude of the k =-12.2 a.u. peak 
obtained here is close to 0.004, and that the same value can be observed in figure [9)d) 
for |* s (fc = ±12,2 a.u.,t = 395 ps)\. 



the photo-association pulse. 
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Appendix: Kinetic energy "gains" in the ground state due to the 
off-resonance cycling of population between the strongly coupled electronic 
states. 

We shall consider the time-dependent Schrodinger equation ([2]) for constant coupling 
Wl- In [28] it is shown that, if the impulsive approximation is valid at a given R, the 
wavefunction on the ground state surface after a time t becomes: 

¥ g A {R,t) =e-^ 9 V^ A '{cos(fit) + icos^sin(fit)}^ g ( J R,0), (A.l) 
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where E g is the energy corresponding to the initial stationnary wavefunction at t=0 
([T + V g ]Vg(R,0) = E g V g (R,Q)), 2A(R) = \V e (R) - V g (R)\ is the local detuning, Q(R) 
the spatial dependent Rabi pulsation: 



Ml(R) = sJWl + A(R) 2 (A.2) 

and 

Ml(fl) Jwl + A(fl)2 



We calculated numerically the time- dependent prediction for the ground state 
wavefunction in the impulsive approximation, ty IA (R, t), using (1A.1I) . In figure lATl we 
show the corresponding momentum distribution at t=300 ps, \^ g A (k, t = 300 ps)\, which 
indeed displays peaks at k values deduced from the maxima of the initial wavefunction, 
according to the discussion of Section 14.31 



In the following we shall introduce approximations into (IA.1I) in order to 
make appear analitically a momentum k associated with a maximum of the initial 
wavefunction at Rq. 

In (IA.ll) . e~^ E9t ^/ g (R, 0) represents a free evolution in the ground state. The factor 
superposed on this free evolution can be separated: 

A(R, t) = e-* At {cos(nt) + i cos 9 sin(fit)} (A.4) 

- r^ 1 + m> + - m> < A - 5 > 

We shall make approximations on the expression flA.5j) in the limit A(R) ^> Wl- Indeed, 
this limit is valid at #0 = 89 a , where A(R ) = hA L /2 « 70 cm" 1 , and W L = 13.17 
cm -1 . Then it is possible to approximate: 



/ W 2 W 2 
hn = A ] jl + — ^A + —. (A.6) 

Consequently, at a given R value where the impulsive approximation is valid, ty g (R,t) 
can be decomposed in two terms: 

w*«(i2,t)+*W(i2,0 (A.7) 
t) = \{l + ^ei&e-fo^R, 0) (A.8) 

¥ ? {R > t] = 47§A e "" [2A+ ^ Vi ^^ ( ^ 0) (A ' 9) 
For A ^> Wl the first term is the dominant one. Indeed, a rough approximation 
TiVL m A implies ^^'(R, t)\ ~ \^ g (R, 0)|, the second term having an incomparable 
smaller contribution ( 4 (^ A |s9a = 0.0087 <C 1). It is the second "small" term which 
offers the explanation for the "k- features" observed in our results. Its Fourier transform 



is: 



q?f) (k,t) = J dRe- ikR ¥ g 2) (R, t) (A. 10) 

= e ~**7 dR A^jA e ^ t+kR ^ a{R ' 0) (A - U) 



Cold molecules formation by shaping with strong light the short-range interaction 24 



with E 2 a = 2A + -7^- ps 2A. We shall take into account only the contribution at the 
integral coming from a small domain of R, Dr, around the point Rq, assuming that 
on this domain *ff g (R,0) is a Gaussian of width 5r centered in i? ; ^g(R,fy\Ro ~ 
ty g (Ro, 0) exp[— (R — i? ) 2 /^]; an d that E 2A (R) depends linearly on R near R : 
E 2A (R) ~ E 2A (R ) + ^§t\ Ro (R - Rq). Then: 

fl rfre-^s^lflo'+^-e-" 2 /^ (A.12) 

where r = R — Rq and the initial domain of integration was taken from i?o — $r to 
-Ro + <5r- F° r a £r sufficiently large, the integral in (IA.12I) becomes a 5 function, such 



as: 



i /I TP 

x27r5(--^\ Ih t + k) (A.13) 

In (IA. 13[) one can see appearing the contribution of an ingoing (fc.Ro = — fc-Ro) plane 
wave e -1 ' 2A b ° t+kR o^ of energy E 2 a(Ro) /% and momentum p = hk [29]. The 5 function 
shows that this momentum p = hk = — \^ t is related to the local energy difference 
between the electronic potentials at Ro, or that dp/dt = — -^"\rq- If one associates 
this momentum to the movement of an ingoing particle of mass fi, p = —fidR/dt, then 
one obtains pj\i = dE 2A /dp, and finally E 2 a = p 2 /2fj,. 

Then indeed (1A.13j) makes appear a momentum k connected to 2A(R ), as we 
observed in our results discussed in Section 14.31 and oriented to small R distances. The 
fact that we first observe an outgoing wave has to be due to the rapid reflection of the 
ingoing waves by the inner wall of the ground state potential. 

An estimation of the wavefunction amplitude \I/^(fc, t)^ for R = 89 ao using 
formula (|A~T3|) gives \^f\k, t)\ Ro=89ao ps 0.001 x 5(\^^\ Ro t + k). This qualitative 
estimation obtained analitically for unprecised t corresponds to the amplitudes 
calculated numerically for k = ± 12.2 a.u., from the dynamics and in the impulsive 
approximation, which are in agreement: \^j:(k = ±12,2 a.u., t = 395 ps)\ ~ \^l A {k = 
-12.2 a.u.,t = 300 ps)\ ps 0.004 (see figure Eb) and figure [AT])- 
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